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Abstract 

We consider the Harper model which describes two dimensional Bloch elec- 
trons in a magnetic field. For irrational flux through the unit-cell the cor- 
responding energy spectrum is known to be a Cantor set with multifractal 
properties. In order to relate the maximal and minimal fractal dimension 
of the spectrum of Harper's equation to the irrational number involved, we 
combine a refined version of the Hofstadter rules with results from semiclas- 
sical analysis and tunneling in phase space. For quadratic irrationals uj with 
continued fraction expansion a; = [0;n] the maximal fractal dimension ex- 
hibits oscillatory behavior as a function of n, which can be explained by the 
structure of the renormalization flow. The asymptotic behavior of the min- 
imal fractal dimension is given by amin ~ const. In n/n. As the generalized 
dimensions can be related to the anomalous diffusion exponents of an initially 
localized wavepacket, our results imply that the time evolution of high order 
moments {r'^), g — > oo is sensible to the parity of n. 
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Originally conceived to describe Bloch electrons in a magnetic field Harper's equation 



+ + 2 cos(27ru;n + z/)0„ = Ecpn (1) 



has become a subject of its own right. The parameter uj is the ratio between the magnetic flux 
per unit cell and the flux quantum, 0„ is related to the wavefunction of the Bloch electron and 
z/ is a phase. Harper-like models are of great interest in a wide variety of physical contexts 
(integer quantum Hall effect 0, superconducting networks 0, electrons in superlattices 0) 
and have given rise to strong interplay between physics and mathematics, as e.g. in the fields 
of semiclassics , non commutative geometry 0, and, quite recently, quantum groups [0. 
Experimental investigations on superconducting networks permit a direct observation of the 
ground state energy as a function of the magnetic flux [Q. Furthermore, the resolution of 
the fine structure of Hofstadter's famous butterfly have been made possible by measuring 
the magnetoresistance oscillations in superlattices 

During the last years there has been great interest in the generalized (Renyi) dimensions 
of the spectral measure of this equation and similar quasiperiodic tight-binding hamilto- 
nians [[TO|-p!^. For the Harper equation it was conjectured that the fractal dimension -Dg=o 
is exactly 1/2 for typical irrational numbers |T5|JT^. Due to the Thouless property P^Pj-pi 



however, there is a simple argument providing strong evidence for Dq being strictly smaller 
than 1/2 [|22|. Recently a statistical theory has been presented to explain the numerically 



observed behavior of Dq as a function of the irrational number uj = [0]n] |Q. Until now 
there is no detailed investigation how the generalized dimensions Dg for a generic q depend 
on the incommensurability uj. However, the generalized dimensions play an important role 
for explaining the anomalous diffusion properties of wavepackets found in incommensurable 
models. Quite recently, the link between Dg and the anomalous diffusion (r^) ~ t^'^'? has 



been made explicit by the relation cxg = -Di-g, which has been verified numerically ^4 . 

In this paper we deal with the problem of the maximal and minimal fractal dimensions 
«min = -D+oo and Omax = -D_oo- This is intended to be a first step to general values of q. 
While Oinin (as Dq) exhibits a monotonic behavior as a function of n, we find that Omax 
shows strongly pronounced even-odd oscillations that can be explained by a refinement of 
the Hofstadter rules [|^. Similar behavior of amin and ttmax has been found for a tight- 
binding Hamiltonian associated with substitution sequences using the trace map approach 



|2|,|26 



The paper is organized as follows: We begin by reviewing the Hofstadter rules (1), then 
we consider the multifractal properties of Harper's equation for the case of the golden and 
silver mean (2) and finally pass to more general quadratic irrational numbers (3), where the 
minimal and maximal fractal dimensions of the related spectrum will be discussed in detail. 



I. HOFSTADTER RULES 

The starting point of our considerations is the hierarchical clustering of the spectrum of 
Harper's equation (|l|). By inspecting the band structure of this equation as a function of u, 
Hofstadter noticed that for a given u it consists of three parts, one central cluster (S) 
and two side clusters (R) that are closely related to the full spectra of Harper's equation 
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with renormalised values of tu, namely R{uj) and S{uj) with 




and 



X I I I < id < h, 

R(uj) = I ^ 2' 

1 It^I if i < < 1 ^ ' 



S{^) = {)V''\ " 0<^<i (3) 

^ ^ ' r 1 ;f 1 < ^ < 1^ V ^ 




where {■} denotes the fractional part (for an illustration see p7| , p!8| 



These clustering rules, denoted as Hostadter rules, have been derived afterwards by 



different methods [IS, 2' 



Until now, however, for practical calculations, it has always been assumed that the side 
clusters and the central cluster are simply rescaled versions of the full spectrum for the 
renormalized uo. Furthermore only the case of the golden and silver mean (o^goid = [0;T], 
f^siiver = [0;2]) has been considered. These particular values of u) are fixed points of both 
maps, R and 5, and therefore self-consistent equations for the fractal dimension Dq can be 
obtained easily [0. 

In order to refine the assumption that the side and central clusters are merely ho- 
mogeneously rescaled versions of the full spectrum for the renormalized a;, we rewrite 
Hofstadter's clustering rules in a more quantitative way. Consider the spectrum of a 
rational approximant uo = p/q with q bands. The partition of these bands on the 
lower side cluster, the central cluster and the upper side cluster is given according to 
den uj = q = den R{ijj) + den S{uj) + den R{lj), where "den " denotes the denominator. 
Therefore we can write the q energy values (upper or lower band edge, for large q, in 
which we are interested, the bandwidths go to zero, so there is no difference) in the following 
form 



'f-{E,{R{u;))) for ^ = 1, ■ ■ ■ , den R{u) 

E,{u) = { PSEi-A.nRi,u.){S{uj))) for I = den R{u) + 1, ■ ■ ■ , den R{uo) + den S{uj) 

/+(Ei_den/?(a;)-den5H(^(t^))) for i = dcu i?(u;) + dcu S'(u;) + 1 , ■ ■ ■ ,g. 

(4) 

The functions f^'^{E) are the fingerprint of Hofstadter's rules and are essential for calcu- 
lating the fractal dimensions of the spectrum for a given uj. The irrational number uj can be 
developped in a continued fraction. Its truncation at the fcth level yields the approximant of 
uo of generation k. An essential point is that f^'^{E) only depend weakly on the generation 
of the approximant of uo. 

For a quasiperiodic tight-binding model where the Hofstadter rules also apply, the func- 
tions f^'^{E) have been calculated in a linear approximation by means of a renormalization 
approach for the case of the golden mean. Since Wgoid is invariant under S and R, the Hof- 
stadter rules close up after the first step, leading to a self-consistent equation for the fractal 
dimensions . 



A generic quadratic irrational number is not invariant under R and S, but successive 
applications of R and/or S involve a finite number of different irrational numbers. This 
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can be seen, if the renormalization equations for uj are rewritten for the continued fraction 
expansion uj = [0; Oi, a2, • ■ ■ ]: 



R{iu) 



[0; a2, as, 04, 
[0; as, 04, 05, 



and 



S{u;) 



'[0;ai -2,02,03, 
[0; 03,04,05, ■■■] 
[0;o2 - 1,03,04, 

^ [0; 04, 05, oe, ■ ■ ■ ] 



if Oi > 1, 
if oi = 1 



if oi > 2, 

if Oi = 2, 

if Oi = 1, 02 > 1, 

if Oi = 1, 02 = 1. 



(5) 



(6) 



Since quadratic numbers have periodic continued fraction expansions, only a finite num- 
ber of different irrationals will occur by applying iteratively R and S to u. Thus, once the 
functions fui{E) are known for all uj involved in the renormalization flow of the initial uj, the 
multifractal properties of the spectrum are determined. 

We will restrict ourselves to quadratic numbers of the form uj = [0;n]. Examples of the 
renormalization flow are given in the following schematic representation for n = 3,4: 



■[0;41: 



R 



R 



[0;1,3] 



[0;2,4] 



R 



R 



[0;2,3] 



R 



The structure of these diagrams depends on the parity of n: For n odd, there are n 
irrational numbers involved and n + 1 cycles occur, n being of the form S^R and one of the 
form S"'. For n even, there are n/2 irrationals and n/2 cycles occur, n/2 — 1 of the form 
S^R and one of the form S""/^. 

Each of these cycles do only leave invariant the exact irrational number uj. Applied to 
an approximant of uj of generation k, an approximant of uj, but of lower generation k' < k is 
obtained, as is immediately obvious by inspecting (|^) and (|^). Thus we can assign to each 
cycle a "generation loss" d := k — k' (cf. Fig. [l|). 
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II. GOLDEN AND SILVER NUMBER 



For the sake of simplicity, we will begin by dealing with the problem for the case 
uj = co'goid,siiver5 which are fixed points of the Hofstadter rules. The mapping f^'^{E) has 
two different fixed points, denoted as R and S in Fig. |], which correspond to the edges re- 
spectively to the center of the spectrum. Assuming in a first approximation that the scaling 
properties of the spectrum are well described by the contraction factors at the centers and at 
the edges, zs = — ^Icenter and zji = ledges, we find that after k steps of renormalization, 
the bandwidhts scale with A ~ z*' while the system size scales with ~ uj~^''^, where uo is 
the algebraically conjugate to uj {uj = ^{n + -\- A) for uj = [0; = \{—n + yjn^ + 4)). 

With A ~ A^~i/" we obtain the two local scaling exponents at the edges and the center 
of the band: 

and 

OiR = , 8 

In zr 

where d{S) and d{R) are the generation losses of the cycles R and S. 

Furthermore, within this approximation a self-consistent equation for r(g) = {q — 1)-Dg, 
the Legendre transformed of /(a), can be given directly: 

2^d{R)q ^d{S)q 

1 ; — = 1, 1^ = ^^gold, t^silver- (9) 

For the case of the golden number, this equation has recently been used to derive the mul- 
tifractal properties of a tight-binding Hamiltonian associated with the Fibonacci sequence 
in the limit of strong modulation |T^. Despite of being a good first approximation, the 
corresponding /(a) has two serious drawbacks: First, /(amin) or /(ttmax) (depending on the 
parameter) does not vanish in contradiction to numerical simulations and indications from 
the trace map approach ||12[. This problem is due to the oversimplification of /^'° by three 
linear pieces, yielding a artificial degeneration of the widths of the bands (2" bands of uj"" 
ones have the same widths), which leads to f{aR) 7^ 0. 

Second, in this approach the extremal scaling properties are found at the band edges or 
at the center of the spectrum. This is, however, by far not the general case. It is true for 
i^goid for both the Harper model and the quasiperiodic Hamiltonian mentioned above, but 
neither for cUgUver; nor for higher irrational numbers = [0; n], n > 3. 

In order to construct an approach eliminating these two problems, we consider the map- 
ping of the spectrum by the renormalization flow as a discrete dynamical system. For a;goid 
and cUsiiver the spectral measure is given by the invariant measure of the map f^'^{E). The 
band edges and the band center are attracting fixed points of f^'^{E), thus giving rise to 
nontrivial scaling behavior as we have seen above. To get accurate approximations of the 
/(a;)-curves, however, we have to take into account higher periodic orbits of the underlying 
dynamical systems. The simplest approximation is to include the orbit of period two. 
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The energies of the two (symmetric) fixed points are given by f^{E) = —E. Therefore we 
find for the corresonding scahng factor 

d{R^R^)\\iuj diR)\\iuj , . 

= -In = (l"' 

where zr^ is the local contraction factor at the fixed points of order two: 

d/±'0(E) 



dE 



(11) 



Taking into account the three contraction factors zr, zs, zr^ (denoted as Sz-model in the 
following) we arrive at the following self- consistent equation for T{q): 



^d{R)q ^d{R)q ^d{S)q 



H 1 — — 1, ^ — l^gold, "^silver- (12) 



^R ^R± 



If zrj^ equals zr equation ([T^) reduces to (H). 

The minimal and maximal dimensions for the 32;-model are given by ayain{ma.x) = 
min(max){ai^, as, By performing the Legendre transformation of (|1^) we find the cor- 

rect behavior /(amin,max) = 0. While for the golden number a^ax = C(s > C(r± > C(r = c^min, 
i.e. the maximal scaling exponent occurs at the center of the spectrum, for the silver number 
the maximal dimension is given by Omax = C(r± > cts > c^r = ctmin, i-e. the bands of maxi- 
mal widths are not at the center of the spectrum, but at the energy given by f^{E) = —E. 
Therefore a linear approximation of f^'^{E) (i.e. two scaling factors, one for the central, 
and one for the side bands) is not sufficent for explaining amin and Omax, even in the case of 

'^silver • 

The approximative /(a)-curves calculated from equation ([12|) are compared in Fig. § 
with those obtained by the usual multifractal formalism |p . Despite of the simplicity of our 
model there is a good agreement between the curves. Including linearization around further 
periodic orbits of the non-linear map f^'^{E) would yield further quantitative improvement. 
We consider the 3z-model as the simplest qualitatively correct approximation. 

Let us note that the Thouless property D^i = 1/2 (i.e. t(— 1) = —1) yields one linear 
equation which the three contraction factors have to fulfill: 

Zr zr^ , _ . _ /,o\ 

H 37^ H 37^ — -1-5 f-^ — "^gold 5 "-^silver {i-'j) 



^d{R) ^d{R) ^d{S) 

Thus, by use of this relation, the calculation of the scaling behavior at the band-edges zr 
and at the center of the band zs would enable to obtain the missing contraction factor zr^ . 



III. HIGHER QUADRATIC IRRATIONAL NUMBERS 

Hitherto we have been considering the case of the fixed points of the Hofstadter rules, 
where the renormalization always yields the same function /^(°)(-E') = f^^^-j{E) = f^'^{E) 
that can be considered as a dynamical system. We have shown that taking into account 
not only the fixed points of this dynamical system, but also a periodic orbit of length two. 



6 



satisfactory /(a)-curves can be obtained. Now we generalize this approach to the irrational 
numbers [0;n] that are not fixed points of the Hofstadter rules but give rise to a nontrivial 
cycle structure that has been shown above. 

In this case we have to deal not only with one function f^'^{E) but with a whole set of 
function {f^f[E)}, where the uji are the irrational numbers occuring in the fiow diagram 
of uj. In order to calculate the scaling properties of the spectrum, we have to consider the 
cycles of the fiow diagram. To each cycle there is one effective f^^{E) that results from 
concatenating the corresponding maps f^f{E) for the Ui occuring in this cycle. For example, 
the effective map of the cycle S'^R in the case of u; = [0; 3] is given by /^^ = 3] °/[o i 3] °'^[o-3] 
(consider equation and the flow diagram). Each of these effective maps can be considered 
as a dynamical system with flxed points and periodic orbits of different lengths. 

Now two natural questions arise: Which cycle of the renormalization flow does lead to 
the minimal (respectively maximal) local scaling exponent? Which flxed point or periodic 
orbit of the dynamical system corresponding to this cycle does yield the extremal scaling 
exponent? In a flrst step, these two questions can be addressed numerically by determinating 
the index of the band of maximal (respectively minimal) width and following its way under 
the renormalization flow. 

We flnd that the the minimal scaling exponent for oj = [0; n] occurs for the cycle R (these 
oj are flxed points of R), i.e. the smallest bands occur at the edges of the spectrum. More 
interestingly, the maximal scaling exponent occurs for the cycle S"'~^R (n odd) respectively 
for S^^^R {n even). The lengths and the structure of these cycles show a striking even-odd 
effect which is a direct result of the Hofstadter rules. Furthermore, for n odd the maximal 
scaling exponent is given by the flxed point of the effective dynamical system, while for even 
values of n the asymmetric orbit of length two yields the maximal scaling exponent. 

For the minimal exponent we flnd for u = [0;n]: 

ttmin — , /r„ —IN ) \^^) 

\nzR{[0;n\) 



and for the maximal exponent we obtain 

21nc:; 



Or, 



In 



/eff'(^r) 



In 



/cff'(^2 



if nodd (15) 
if n even, (16) 



where E^ {E2) is the energy of the flxed point (energy of the orbit of period two) of the 
effective map f^{E). Replacing the effective map f^{E) by the individual maps, we flnd 

2 IncD Til /-.rrX 

Ei=M^2 In ^5 ( [0 ; 2 , n] ) + In ( [0 ; 2 , n] ) 
amax = -;™ i --nil nn o -n " neven. (18) 

Ei=4,i even Zs{[0; I, Tl]) + lu Zr^ ( [0; 2, n] ) 

Inspecting these two equations one is lead to the assumption that Omax as a function of 
u = [0;n] should be sensible to the parity of n, while ttmin is expected to show a monotonic 
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behavior function of n. This is indeed the shows Fig. ^. For a more quantitative 

discussion of our approximation, however, the behavior of zr, zr^ and zs as a function of uj 
is required. For hmiting cases, as — >■ or o; — ^, the contraction factor zs can be obtained 
by semi classical analysis and by considering tunneling in phase space [^ , ^ ,p|,plH33[| , if one 



makes the further assumption that f^{E) is nearly linear, so that zs can be written as ratio 
of band edges: zs{uj) = Ei_{uj) / Eq{S (u)) . For u ^ zs{uj) is therefore given by the ratio 
of two Landau levels and with the results of fHl we find: 



zs{uj) = l-7iu + nu^ + (27r - vr^ + — )cu^ + 0{uj^) {u 0). (19) 
For a; — > i tunneling has to be taken into account for determining Ei_{lj) and we find 

i bsil-uj) + cs{l-ujf + ds{l-u)Hn{l-iu) {uj^l), (20) 



Inzsiu) '2 ' '2 ' '2 ' '2 ' ' 2' 

where 65 is given analytically by 65 = vr • (/J^^^ ln(cosA: + Vl + cos^ k)dk)-^ ^ 3.429815 and 
cs and are fit parameters (cf. appendix). Although equations ([T9|) and ( pOf ) are only 
expansions near u = and u = ^, they describe ^;5(ci;) over the whole range < to < ^ 
amazingly well, if one takes equation ([19|) for < < 0.25 and equation ([20| ) for 0.25 < 
u; < 0.5 (cf. Fig. |). 

In principle, a similar analysis can be performed for the mean slope = {Eo{uj) — 
Ei^{uj)) / {2Eq{R{uj))) . As the linear approximation for f^{E) assumed in this relation is a 
rather crude one for general uj, the values of zr and zr^ can deviate considerably from Yr. 
This, however, is not the case for — >• 0, where zr and zr^ are well approximated by 
that is given in this regime by 

- , , = hRUj + cruj'^ + d.Ru'^ In cu, (tu 0) (21) 

\nzRiuj) 



where 6^ = 27r ■ (/g ln(2 - cos A; + v/(2-cosA;)2 - l)dA;)-i 0.857454 and cr ^ 1.5 and 
(Ir ^ 0.37 again are fit parameters. 

As amin is determined by zr[0; n] ^ ^^[0; n] (as [0; rl] — > forn — 00), we insert equation 
(pl|) into (|T4|) and obtain 

ctmin — {bRU + Cruj"^ + dRUj"^ \tiuj) Incj, (22) 

where uj = |(v^^M-4 — ra) and uj = ^{n + \/rfi + 4). Therefore, asymptotically amin ~ 
feijlnn/ra (ra ^ 00) and lim„^oo ttmin = 0. 

For Omax the values zr and z/?^ for uj near to ^ are needed. As they are not directly 
available in our approach, we take approximately zr{iJ) = TrIuj) and zr^ = 1.4^^(0;), where 
Yr^uj) is given by equation (pi]). This approximation is intended to give the qualitative 
behavior for Omax- 

In Fig. ^ the predictions for and Omax for 3 < n < 10 are compared with the corre- 
sponding values obtained by diagonalization. For a^nin and n > 5 there is good agreement 
between pediction and simulation, while for n < 5 the predicted value of amin are to high, 
reflecting the fact that zr < Tr. For ctmax we find qualitatively the even-odd oscillations, but 
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of course, due to the problems mentionned above, we can not hope to achieve quantitative 
agreement. In particular an open question is the behavior of Omax for n — ^ oo. For the 
moment, we can only state, that if In 2:^(0;) and Ihzr^^uj) is bounded for u ^ ^, then ctmax 
does not reach the value 1 for n — > 00. 

Summing up, we have used a refined version of Hofstadter rules together with results 
from semiclassical analysis to link the maximal and minimal fractal dimension of the spec- 
trum of the Harper equation to the continued-fraction expansion of the irrational number 
characterizing the flux per unit-cell. We have shown that the oscillatory behavior of Omax 
for uj = [0;n] as a function of n can be qualitatively explained by the discrete flow diagrams 
of the Hofstadter rules. For an explicit expression has been given which describes fairly 
well the data obtained by diagonalization. We conjecture that the result for Omax can be 
generalized to quadratic numbers of the continued fraction expansion lu = [0; rii, n2, ■ ■ ■ , Up]. 
Take for example u = [0; nin2]- In this case the renormlization flow diagramm contains 
|(ni + different irrational numbers, if both rii and n2 are even. In all the other cases, 
there are rii +n2 different irrational numbers. Thus for even values of ni the cycle structure 
depends on the parity of ^2, while this is not the case for rii odd. Hence, for fixed even 
values of ni we expect even-odd oscillations of Omax as function of n2, while for odd values 
of ni no such oscillations are expected. 

Finally we note that the applicability of the Hofstadter rules to a wide class of quasiperi- 
odic models suggests that our method could be useful for other models than the one we have 
considered here. 



IV. APPENDIX 

This appendix is intended to provide the calculations for zs and in the semiclassic 
approximation using methods and results described in [^,p|,plH33[ . 



Equation ([T9| ) can be obtained directly by inserting the semiclassical expansion of the 



Landau levels up to order three [§] into 



For UJ ^ the contraction factor Zs is given by 

^0(4(2 -a;)) 2 

While the denominator corresponds to Landau levels near u = and the already mentioned 
semiclassical expansion can be used, the numerator involves the calculation of i?i_(a;) for 
u ^. Since for = | the effective Hamiltonian H = ±2 V cos^ ki + cos^ k2 has the form 
of a Dirac operator at the four critical points (±f , ±f ) (independent signs) we are left 
with the calculation of tunneling between Dirac levels. 

Although the usual formula for level splitting due to tunneling 

AE ~ const -h-exp (-^—Y) ' (^^^ 
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where S is the action integral between degenerate minima, h is Planck's constant and the 
proportional constant is related to the classical frequency, was originally derived for tunneling 
between parabolic minima |Q, we tentatively use it for our case, where the minima are 
conical. With h = 2tt{^ — u) we obtain 

B._M~(l/2-.)exp(-^^^f^^). (26) 

where the action integral has to be calculated between two adjacent minima in phase space. 
The (complex) orbits are given hj H = 0, i.e. cos^ ki + cos^ k2 = and therefore 

S= arccos(±i COS A;i)dA;i (27) 

J-n/2 

|Im^| = 2 / ln(cos fci + V 1 + cos^ A;i)dfci ^ 1.831931. (28) 



Formula (|2q) with (pSD inserted is compared to numerical results from diagonalisation in 
Fig. 1^. We note that the fact of the action being purely imaginary shows the absense of 
braiding of the Dirac levels. To calculate zs[io — 1/2) we insert (^6]) together with 

£;o(4(l/2 - cj)) ~ 4 - 87r(l/2 - cj) + O ((1/2 - ujf) (29) 

into (p^) and obtain (pO|). The parameters cs ~ 9.6 and ds ~ 5.2 are determined by fitting 
to the numerical data. 

The calculations for T^iio 0) are analogous, but only Landau levels have to be con- 
sidered. 
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FIG. 1. Cycles occuring in the discrete renormalization diagram for lj = [0;n]. 
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Eo(S(co)) 

E 



Ei_(co) Ei^(co) Eo(a)) 

f--°(E) 

FIG. 2. Schematic representation of f^'^{E). The points R and S are the fixed points of 
f^'^{E). The point R± corresponds to the two-orbit E — > —E —>■ E. 
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a 

FIG. 4. The multifractal spectrum /(a) for the Harper model for Wgoid (narrow curve) and 
"^silver (large curve) compared to the 3z-model. The /(a)-curves obtained from diagonaliza- 
tion and the thermodynamic formalism are drawn with symbols, the f{a) obtained from Leg- 
endre- Transformation of (|l2|) are depicted with lines (zf'*^ = 0.103, zf^ = 0.164, zf^"^ = 0.0718; 
^silver ^ 0.118, zf2' 



^silver ^ 0.118, zf^^' = 0.248, 4"™'' = 0.020). 
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FIG. 5. The contraction factor 2:5(1^) = Ei^{uj)/Eq{uj) compared to the results from semiclas- 
sical analysis (equation (19) for < lo < 0.25 and equation (pO| ) for 0.25 < lj < 0.5 (65 = 3.429815, 
cs = 9.6, ds = 5.2)). Note that zs{io) is remarkably smooth compared to Ei^{u)) and Eo{u) (Fig. 
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FIG. 7. The broadening of the Dirac level at = 0, a; = ^ compared to the result from 
tunneling in phase space — l/ln£?i_(a;) ~ 27r(^ — a;)/|ImS'|, where the imaginary part of the action 
is given by [ImSI = /^^^2 ^^(^°^ A' + Vl + cos^ k)dk ^ 1.831931. For 1/2-u < 0.013 the broadening 
of Dirac levels can not be resolved numerically. 
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